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The density functional theory of nuclear structure provides a many-particle wave function that is 
useful for static properties, but an extension of the theory is necessary to describe correlation effects 
or other dynamic properties. Here we propose a procedure to extend the theory by mapping the 
properties of the self-consistent mean-field Hamiltonian onto an effective shell-model Hamiltonian 
with two-body interactions. In this initial study, we consider the sd-shell nuclei ^°Ne, ^^Mg, ^®Si, and 
Ar. Our first application is in the framework of the USD shell-model Hamiltonian, using its mean- 
field approximation to construct an effective Hamiltonian and partially recover correlation effects. 
We find that more than half of the correlation energy is due to the quadrupole interaction. We then 
follow a similar procedure but using the SLy4 Skyrme energy functional as our starting point and 
truncating the space to the spherical sd shell. The constructed shell-model Hamiltonian is found to 
satisfy minimal consistency requirements to reproduce the properties of the mean-field solution. The 
quadrupolar correlation energies computed with the mapped Hamiltonian are reasonable compared 
with those computed by other methods. 


PACS numbers: 21.60.-n, 21.60.Jz, 21.60.Cs, 21.10.Dr 

I. INTRODUCTION 

Quantitative theory of nuclear structure has followed 
two different paths. One path is a self-consistent mean- 
field theory (SCMF) 0, also called density functional 
theory,^ and the other is the configuration-interaction 
shell model (CISM) approach Q. The former requires 
only a few parameters and has broad applicability across 
the table of nuclides, while the latter is specific to re¬ 
gions of the nuclides and as usually formulated requires 
many parameters to specify the Hamiltonian. However, 
the SCMF is of limited accuracy for calculating binding 
energies, and the treatment of excitations requires exten¬ 
sions of the theory. In contrast, the CISM gives a rather 
accurate and almost^ complete description of the low- 
energy spectroscopy. Brown and Richter ^ proposed to 
use the SCMF to determine the single-particle Hamilto¬ 
nian of the CISM so as to gain advantages of both ap¬ 
proaches. In this work we go farther, using the SCMF to 
also determine parts of the two-body interaction Hamil¬ 
tonian. 

Our work has several motivations. One is to develop 
a well-justified and accurate theory of nuclear masses. 
SCMF is systematic with a strong theoretical foundation, 
but it may have reached its limits of accuracy at a level 
of an MeV rms residual error over the table of known 
masses. Some effects of correlations are missing from the 


^ This terminology is used because of the similarity to the density 
functional theory of many-electron systems Q. 

^ The qualification is due to the common occurrence of intruder 
states. 


theory, and indeed better accuracy can be achieved when 
phenomenological terms are included to take them into 
account ■ A systematic calculation of those correlation 
energies would be possible if we could find a way to map 
the SCMF on a shell-model Hamiltonian^. 

Another motivation is to bring more precision to 
shell-model calculations of one-body observables such 
as quadrupole moments and transition intensities. For 
example, the SCMF is well adapted to computing 
quadrupole matrix elements because of the large single¬ 
particle space that can be treated. The CISM requires 
severe truncation of the space, and a systematic mapping 
procedure would determine the effective operators to be 
used in the smaller truncated shell-model space. 

Finally, such a mapping would be helpful to construct 
a global theory of nuclear level densities. While the over¬ 
all behavior of level densities can be derived from the 
independent particle shell model, there are significant in¬ 
teraction effects. The shell model Monte Carlo (SMMC) 
method has been efficient for calculating level densities in 
large shell-model spaces |E S 0 > but it requires as input 
a parameterized one- and two-body shell-model Hamil¬ 
tonian. It should be mentioned that the SMMC requires 
for its efficiency interactions that have a specific sign.^ 


® The SCMF energy functionals in use are fitted to the binding 
energies, and in principle they should be refitted if the correlation 
energies is added as a separate contribution. In practice, the 
refit does not change the parameters significantly so the original 
functional may be used in the determination of the Fock-space 
Hamiltonian. 

Interactions containing small terms havi^ the bad sign can be 
treated as well with the method of Ref. . 
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These conditions are satisfied for the Hamiltonian we will 
construct here. 

Here we propose a mapping procedure and apply it to 
sd shell nuclei. This region of the table of nuclides is a 
paradi gm for the CISM. For example, the USD Hamil¬ 
tonian defined in the sd shell space reproduces very 
well binding energies, excitation energies, and transition 
rates between states. For the actual SCMF that one 
might wish to map, we will use one of the Skyrme family 
of energy functionals. The corresponding SCMF equa¬ 
tions can be solved with computer codes that are publicly 
available PIH; in particular we shall use in this work 
the Brussels-Paris code mm, together with the SLyd 
parameterization of the Skyrme energy functional |l4l| . 
This functional is fit to a number of empirical properties 
of nuclei and with a slight readjustment gives a 1.7 MeV 
rms fit to masses of ~ 600 even-even nuclei [Hill. 

The outline of this paper is as follows. In Section QH 
we discuss the main idea and procedure for construct¬ 
ing the SCMF to CISM map. Our present study fo¬ 
cuses on quadrupolar correlations. In Section mil we test 
our mapping procedure using the mean-field approxima¬ 
tion of the USD Hamiltonian as the SCMF theory to be 
mapped on a shell-model Hamiltonian. In Section ESI we 
discuss the technical construction of the map, starting 
from the SCMF theory of the Skyrme family of energy 
functionals (in particular SLy4). In Section 0 we carry 
out the mapping explicitly for the nucleus ^®Si, while in 
Section rVn we summarize similar results for all other de¬ 
formed N = Z sd shell nuclei. In Section we use 
the mass quadrupole operator to construct an effective 
quadrupole-quadrupole interaction. Finally, in Section 
Enn we discuss possible future developments of the the¬ 
ory introduced here. 

II. METHODOLOGY 

Our goal is to choose a suitable truncated shell-model 
space and construct within it a spherical shell-model 
Hamiltonian containing one- and two-body terms. We 
write this CISM Hamiltonian in the form^ 

H = ~ ^ ■ ( 1 ) 

cx.^j^m K 

Here Caj are single-particle energies of spherical single¬ 
particle orbitals. The orbitals are characterized by the 
quantum numbers ajm, where j,m is the total single¬ 
particle angular momentum and its projection, and by 


® Throughout this work we will use a consistent notation with 
the following conventions. Operators in Fock space (i.e., opera¬ 
tors containing particle creation and annihilation operators) are 
distinguished with a circumflex. Matrices representing single¬ 
particle operators in a single-particle space are denoted in bold¬ 
face. The boldface is dropped for individual elements of those 
matrices. 


a, denoting the remaining single-particle quantum num¬ 
bers®. Ai^ are one-body tensor operators and are in¬ 
teraction strengths. The index k incorporates all the 
quantum numbers necessary to specify a spherical tensor 
operator of definite rank K^. In general 

A^m = Ao'i' , 

ajcx'j' mm' 

( 2 ) 

and the scalar product^ in Eq. o is Ai^ ■ Ak, = 

Eq. O provides the most natural form to which the 
SMMC method can be applied. With enough terms in 
the K sum, it is sufficiently general to treat any two- 
body interaction that is rotationally invariant and time- 
reversal symmetric. The maximal number of terms k for 
any given rank of the multipole operators is the 

number of orbital pairs whose angular momenta j^j' 
can be coupled to K^, by angular momentum and par¬ 
ity selection rules. We also note that faj,a'j' 3''^® P''0" 
portional to the reduced matrix elements of i.e., 
faj,a'r = {oijWAAWj')/+ 1. In the development 
below we do not need to distinguish orbitals with differ¬ 
ent quantum numbers a, and we will drop that index. 

Our starting point is the SCMF, which minimizes an 
energy functional for a configuration of single-particle or¬ 
bitals (j)k (i.e., a Slater determinant) or the BCS general¬ 
ization thereof. Besides the ground-state minimum, one 
can obtain constrained solutions by minimizing the en¬ 
ergy functional in the presence of an external field XQ, 
where Q is the mass quadrupole operator defined below 
and A is its strength. The SCMF delivers a set of orbitals 
(j)^ and corresponding single-particle energies e^, as well 
as the total energy E^. Our task is to use this infor¬ 
mation to construct a Hamiltonian of the form (^. The 
minimal requirement of such an effective Hamiltonian is 
that it reproduces ground-state SCMF one-body observ¬ 
ables sufficiently well within a mean-field approximation 
(in the truncated space). We use this requirement as a 
guide to construct the effective Hamiltonian. 

In the applications discussed in this work, we treat 
nuclei whose SCMF ground state is deformed. This will 
simplify the mapping to some extent because we do not 
have to impose a constraining field. 

The first step is the construction of a spherical single¬ 
particle basis. Due to the rotational invariance of the the 
energy functional, a spherically symmetric density matrix 
at a given iteration step leads to a one-body Hamiltonian 
that is also spherically symmetric. For a closed sub-shell 


® E.g., a = tz,n,l with tz the isospin (proton or neutron), n the 
radial quantum number and I the orbital angular momentum 
^ The tensor operator Ak. satisfies with c = 

±1 and is either hermitian (c = 1) or anti-hermitian (c = —1). 
Therefore, the corresponding interaction term in m can also be 
written as Af, ■ Af, = cJ2m A^^^Ai^m- 
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nucleus (e.g., ^®Si) the density matrix in the next iter¬ 
ation step will remain spherically symmetric. Thus, for 
a closed sub-shell nucleus, if we start from a rotation- 
ally invariant single-particle density and iterate, we will 
eventually converge to a spherical solution. For an open 
shell nucleus, we have to choose which m orbitals of the 
valence partially occupied j-shell are filled. Such a choice 
will break the spherical symmetry of the density and lead 
to a deformed solution. To overcome this problem, one 
may use the uniform filling approximation in which the 
last n particles are distributed uniformly over the cor¬ 
responding 2j -I- 1 orbitals (i.e., we assume occupations 
of = n/{2j -I- 1)). As long as this approximation is 
used in each HF iteration, we will converge to a spherical 
solution. In practice, a pairing interaction or BCS-like 
occupations with a fixed gap are used to guarantee the 
spherical symmetry of the solution. 

We will discuss the construction of a spherical basis in 
more detail in Sect ion ITVl but for now let us assume that 
we have found the spherical solution with orbitals (pjm as 
well as the (deformed) ground-state solution with orbitals 
(j)k- The single-particle density matrix of the deformed 
mean-field solution is given by p^f with Pmf,fefc' = nkSkk ', 
and rife is the occupation number of orbital k in the 
ground-state basis. In the HF approximation rife = 1 
for the lowest N orbitals and rife = 0 for all other or¬ 
bitals (here and in the following we use N to denote the 
number of particles of a given type, i.e., either protons 
or neutrons). We transform single-particle operators to 
the spherical basis using the transformation matrix U 
defined by taking the overlaps, Ukjm = Thus 

the single-particle density matrix p in the spherical basis 
is given by 

Pjm,j'm' — ^ • (3) 

k 

As long as the spherical basis (j)jm is complete, U is uni¬ 
tary and is just a different representation of the same 
ground-state density matrix p^f. However, in construct¬ 
ing the mapping we choose a truncated model space, and 
the matrix U (with jm values restricted to this space) is 
no longer unitary. This p is then the density matrix in the 
truncated space, i.e., Pp^fP, where P is the projector on 
the truncated single-particle space. We can calculate ex¬ 
pectation values of one-body observables in the truncated 
space by using this truncated density Ppj^^fP. For exam¬ 
ple, the expectation value of the mass quadrupole op¬ 
erator Q = in the truncated 

space may be computed as 

{Q) = tr (qPpjjjfP) = ^ ^ qjmj'm' PjmJ'm' > (4) 

where qjmj'm' = {jiTi\Q\j'm') are the matrix elements of 
the mass quadrupole operator Q = — xj — yf) in 


the single-particle space.® 

Since the transformation U from the deformed basis 
into the truncated spherical basis may not be unitary, 
it is important to choose the truncated single-particle 
space such that the truncated ground-state density ma¬ 
trix PpjjjfP satisfies some minimal properties required 
for a density matrix that is the solution of the CISM 
Hamiltonian in a mean-field approximation. One such 
requirement is that the number of particles be correct; 
this is checked with the formula 

tr p = TV . (5) 

This property is satisfied exactly in the full space with N 
being the total number of particles, and the truncated p 
should satisfy this condition approximately with N being 
the number of valence particles. Another requirement is 
that the density matrix represent a single configuration 
(i.e., a Slater determinant). This can be expressed in the 
matrix equation 

pp = p . (6) 

Equivalently, the requirement is that the eigenvalues of 
the density matrix be zero or one (with exactly N eigen¬ 
values equal to one). We will examine the eigenvalue 
spectrum of our mapped density matrices to see how well 
this requirement is satisfied. 

The mean-field Hamiltonian in the SCMF ground- 
state basis is the single-particle matrix hmf with elements 
hmf,kk' = efc<5fe,fe'. Its matrix representation in the trun¬ 
cated spherical basis is given by 

h = Uth^fU . (7) 

Let us expand h in multipoles, writing h = 
where is an irreducible tensor of rank K. The re¬ 
duced matrix elements of are given by 

^ m,m' 

( 8 ) 

We now assume that the SCMF ground state is axi¬ 
ally symmetry, so that only the M = 0 components of 
h(^) are non-vanishing. We shall also restrict ourselves 
to even-even nuclei, which by time-reversal invariance 
permits only even multipoles K. The lowest multipole 
is the monopole from which we obtain the spherical 
single-particle energies 

"‘jmj’m’ CjfUrnm' ■ 

These single-particle energies are given by e® = 

Sm ^jm,jm/ (2j + 1). 


In the following we will omit the projector P in the definition of 
the truncated p and it should be understood from the context 
whether p denotes the complete or projected density matrix. 
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When the ground state is deformed, the K = 2 part 
of the single-particle Hamiltonian will be substantial and 
can be used to determine an effective interaction in that 
quadrupolar channel. Defining the second-quantized one- 
body tensor operators (see Eq. O) 




JJ 


^ v2j + 


( 10 ) 


we consider the following effective CISM Hamiltonian 


H = . ( 11 ) 

The Hamiltonian dJ has the form o and m with a 
single multipole term K^. = 2 and f = . 

Assuming axial symmetry, the Hartree Hamiltonian of 

(tTTll is given by hn = — g(hg^^)hQ^^ (where (hg^^) = 

( 2 ) 

tr (phg ')). Choosing g to satisfy 

^ , (12) 


the Hartree Hamiltonian is hjj = -|- hg^^ If this 

procedure is carried out in the complete space (i.e., no 
truncation), then the Hartree mean-field Hamiltonian of 
(El approximately coincides with the mean-field Hamil¬ 
tonian of the original SCMF theory, and their ground- 
state density matrices coincide. If g is chosen to satisfy 
El in the truncated space, then the Hartree Hamilto¬ 
nian of El in the truncated space is approximately the 
SCMF mean-field Hamiltonian projected onto the trun¬ 
cated space. This will be strictly correct to the extent 
that calculated with (i.e., the SCMF den- 

sity in the truncated space) coincides with (hg ') cal¬ 
culated with the self-consistent density matrix p of the 
Hartree problem in the truncated space. 

Another important quantity is the total energy of the 
system. The absolute SCMF energies need not be re¬ 
produced by the mapped Hamiltonian, but the energy 
difference between the mean-field minimum E^{ and the 
spherical state Fgph should come out the same. We call 
this the deformation energy E^ef, 


Edef = Amf — Esph- (13) 

A deformation energy can also be calculated by solving 
the CISM Hamiltonian El in the mean-field approxima¬ 
tion. The above procedure to determine g does not guar¬ 
antee that the deformation energy of the original SCMF 
theory is reproduced by the mapped CISM Hamiltonian. 
We fine tune the value of g so as to match the deformation 
energies of the SCMF theory and of the CISM Hamilto¬ 
nian. As shown below, the deformation energy of the 
CISM Hamiltonian is very sensitive to g, and a relatively 
small change from its value determined by El will suf¬ 
fice to reproduce the correct deformation energy. On the 
other hand, the quadrupole moment and the ground den¬ 
sity matrix (in the truncated space) are less sensitive to 
g and are still expected to match approximately. 


We have used the Hartree approximation to motivate 
our choice of the effective interaction El- However, in 
practice we include the Fock term in the Hamiltonian 
to determine the density matrix and deformation energy 
of the system. Using a = {jm) to denote the spherical 
orbitals, the Hartree-Fock Hamiltonian of El is given 
by 

Kc = h^aJ-gEMi-rh^Mlc {EtAMPdb) 

E9j2M(~)^'l2bd^Ai-ad^-M-bcPdb^ ( 14 ) 
( 2 ') ^ ( 2 ^ 

where = {a I ^A/ 1^)' the menn field solution Pac 

in El possesses axial symmetry, only the M = 0 com¬ 
ponent of contributes to the direct term, while all 
values of M can contribute to the sum in the exchange 
term. 

There are several sources for the imperfections in the 
mapping. Besides the truncations in multipolarity and 
shell orbital space, the underlying SCMF typically has 
three-body interactions that are ignored in the mapped 
Hamiltonian. Thus, it is important to check the mutual 
consistency of the mean-field observables if the mapped 
Hamiltonian is to be used with confidence. 


III. A WARM-UP EXERCISE: THE USD 
HAMILTONIAN 

As a warm-up exercise for the mapping onto the Hamil¬ 
tonian Eq. El, let us consider the USD Hamiltonian El 
as the target and its mean-field approximation as the 
energy functional to be mapped back onto an effective 
CISM Hamiltonian. As mentioned earlier, the USD 
Hamiltonian gives a very good description of the proper¬ 
ties of sd-shell nuclei, including binding energies, spec¬ 
tra and transition moments. We describe our proce¬ 
dure again but now taking the nucleus ^®Si as an ex¬ 
ample. First, we show the energy landscape for the so¬ 
lutions of the constrained HF equation in Fig. ^ One 
sees the ground-state minimum at an oblate deformation 
(Q) « —70 fm^ and a stationary point at zero defor¬ 
mation. In fact we only need the properties of the HF 
solutions at these two points. 

The spherical solution provides us with the spherical 
basis states for expressing the deformed orbitals. In 
the sd shell-model space this is trivial; the spherical 
single-particle HF states coincide with the original single¬ 
particle states of the USD Hamiltonian and only their en¬ 
ergies are shifted.® In particular, a large monopole con¬ 
tribution to the single-particle energies is automatically 


® The spherical state is always a stationary solution of the mean- 
field equations (in the uniform filling approximation), but it may 
not be lowest state satisfying the constraint (Q) = 0. In that 
case, one has to find it in some other way. In this section we have 
found the spherical solution by starting from a spherical density 
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subtracted when the spherical HF energies are calculated. 
The spherical state has a HF energy of Egph = —126.03 
MeV and zero quadrupole moment. 



constrained Hartree-Fock solution of the USD Hamiltonian. 

The deformed mean-field ground state is obtained by 
iterating the HF equations starting from a deformed 
oblate state and without a constraining field. The so¬ 
lution has an energy of E^f = —129.61 MeV and a 
quadrupole moment of {Q) = —20.3 6^ = —69.6 fm^ 
(here b is the oscillator radius).The deformation en¬ 
ergy is the energy difference between the two solutions, 
Edef = 129.61 - 126.03 = 3.58 MeV. The transfor¬ 
mation matrix U is just the matrix describing the de¬ 
formed mean-field orbitals in the shell representation. 
Since there is no orbital space truncation in the mapping 
sd —>■ sd, the transformation matrix U is unitary and 
the consistency requirements on the ground-state density 
matrix © and m are automatically satisfied. 

The next step is to express the single-particle Hamilto¬ 
nian of the deformed solution in the shell representation, 
Eq. 0 , and carry out the multipole decomposition keep¬ 
ing only the K = 0 and K = 2 terms. This determines 
the effective Hamiltonian Eq. m except for an overall 
strength parameter g. 

A simple estimate of g can be obtained from the naive 
Hartree theory, g = = 0.0218 MeV/5^ (see 

Eq. (I12ll b We tune g to match the deformation energy 
of the USD Hamiltonian. In practice we solve the HF 
equations for the effective CISM Hamiltonian Eq. H1 1 1I for 
a range of values of g, making tables of the deformation 
energy E^efig) and the quadrupole moment (Q). We 
then determine a value for g that fits the deformation 
energy of the USD Hamiltonian E^ef = 3.58 MeV. We 


and using the uniform filling approximation at each iteration 
step. 

In the conversion from b to fm we have used = himujo with 
Jtujo = 45 ^- V3 _ 25A-2/3. 


find g = 0.0232 MeV/6"^ and the deformed state has a 
quadrupole moment of -18.2 6^, about 10% lower than 
the actual minimum of the full USD Hamiltonian. 

TABLE I: SCMF with the USD Hamiltonian and the mapped 
Hamiltonian Eq. o for 4A nuclei in the sd shell. Energies 
are in MeV, (Q) is in {b being the oscillator radius) and g 
is in MeV/b^ 


Nucleus 

Interaction g 

Esph -^mf -E^def {Q} Ecorr 

2°Ne 

USD 

h(2) 0.051 

-25.46 -36.38 10.92 15.4 4.1 

-31.42 -42.34 fit 15.0 3.5 

24 Mg 

USD 

0.0261 

-68.20 -80.17 11.97 18.0 6.9 

-97.04 -108.98 fit 17.6 4.4 

to 

00 

W 

USD 

h(2) 0.0232 

-126.03 -129.61 3.58 -20.3 6.3 

-182.11 -185.72 fit -18.2 4.3 

36 Ar 

USD 

h(2) 0.062 

-222.75 -226.56 3.82 -13.5 4.0 

-372.75 -376.55 fit -12.5 1.9 


We have made similar calculations for the other de¬ 
formed N = Z even-even nuclei in the sd shell. This 
excludes ^^S, which is found to be spherical in the HF 
approximation. The results are summarized in Table 
I. We determine g by matching the deformation en¬ 
ergy of the USD Hamiltonian, and in all cases we find 
the quadrupole moment of the CISM Hamiltonian to be 
within ^ 3% — 10% of the quadrupole moment of the 
USD Hamiltonian (where both moments are calculated 
in the HF approximation). 

In Fig. 13 we show the occupations {fij) = of 

the spherical orbitals versus the spherical single-particle 
energy in the deformed AN sd-shell nuclei. The occu¬ 
pation numbers for the effective Hamiltonians are com¬ 
parable to the occupations for the USD interaction. 

An important use of the effective CISM Hamiltonian 
will be to calculate the correlation energy Ecorr, defined 
as the difference between the mean-field energy E^{ and 
the CISM ground-state energy Egs, 

Ecorr — Eynf ^gs • (1^) 

The quadrupolar correlation energies extracted from the 
effective CISM Hamiltonians are plotted versus mass 
number A in Fig. 13 (solid squares) and compared with 
the full USD correlation energies (open circles). Note 
that the quadrupole-quadrupole interaction is responsi¬ 
ble for more than one half of the total correlation energy 
in all cases but ^®Ar. 


IV. MAPPING THE SCMF 

The SCMF of interest are based on local or nearly lo¬ 
cal energy functionals and use large bases to describe the 
single-particle orbitals. This introduces technical prob¬ 
lems in determining the spherical orbitals as well as con¬ 
ceptual problems associated with the truncation of the 
space. This section is devoted to these issues. Some of 
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FIG. 2: Occupation numbers (fij) of the spherical orbitals in 
the deformed mean-held solution as a function of the spherical 
single-particle energy in four sd-shell nuclei (^°Ne,^'*Mg, 
^®Si and ®®Ar). Results for the effective Hamiltonian I1 1 1I 
(solid squares) are compared with results for the USD inter¬ 
action (open circles). Also shown are results for the effective 
interaction discussed in Section Em (solid triangles). 



FIG. 3: Correlation energies iScorr versus mass number A in 
deformed N — Z sd shell nuclei. The symbols are as in Fig.|^ 
Quadrupolar correlations are responsible for more than 50% 
of the full correlation energies (with the exception of ®®Ar). 


the details are specific to the ev8 SCMF code that 
we use for our numerical computations. 

Since the effective CISM Hamiltonian H must be ex¬ 
pressed in a spherical shell basis, we need to find the 
spherical SCMF solution to construct the basis. We will 
denote the single-particle wave functions of this spherical 
solution by ipajm and truncate the space according to the 
needs of the CISM. Whether or not a solution is spher¬ 
ical can be easily determined by the 2j + 1 degeneracy 


in the single-particle energies. The (Q) = 0 HF ground 
state may or may not be spherical. Usually one includes 
a pairing field in the SCMF, which stabilizes a spheri¬ 
cal solution for (Q) = 0. Our starting wave functions 
are obtained from the archive M , which were calculated 
including a delta-function pairing interaction with a re¬ 
alistic strength. We take the SCMF wave function corre¬ 
sponding to {Q) = 0 (which is spherical), and perform a 
few unconstrained iterations with the pairing interaction 
turned off. However, we include BCS-like orbital occupa¬ 
tion factors vf = {1 — {ck — m)/\/( ffc ~ fJ-)^ + ^^)/2 with 
a fixed pairing gap A (/r is the chemical potential). This 
procedure is applied for several values of A, e.g., 2, 1 and 
0.5 MeV. The resulting total spherical energies are found 
to have an approximate linear dependence on A and we 
have used a linear extrapolation to estimate the value at 
A = 0. We take this extrapolated value to be the HF 
energy of the spherical solution. 

We use the code ev8, which assumes that the SCMF 
solution has time-reversal symmetry. The orbitals then 
come in degenerate time-reversed pairs, so only half of 
them are considered, e.g., the positive z-signature^^ or¬ 
bitals [T^ . For an axially symmetric solution, Jz = m 
is a good quantum number and the positive z-signature 
orbitals have m = 1/2,—3/2, 5/2,—7/2,_ As a prac¬ 

tical matter, there may not always be a clean separa¬ 
tion between different m states (belonging to a given j- 
shell) and even between orbitals of different j-shells, due 
to degeneracies and numerical imprecision. In practice, 
we first determine j by identifying, within the manifold 
of positive z-signature states, the j -b 1/2 nearly degen¬ 
erate states with a given parity. We then diagonalize 
the (j + 1/2) X (j + 1/2) matrix representing the axial 
quadrupole operator Qq = Q to obtain states of good 
\m\. Using the Wigner-Eckart theorem, we have 

{jm\Q^Q^\jm) = Cj{j\\Q^'^'>\\j)[3m^ - j{j + 1)] (16) 

with Cj = 2[(2j — l)2j(2j -b 1)(2/ -b 3)]“^/^. Thus by 
arranging the eigenvalues of Q in ascending order, we 
identify the orbitals with \m\ = j, j — 1 ,..., 1/2.^^ 

Each good-|m| orbital constructed above is determined 
up to an overall phase. However, in the spherical shell 
model, the relative phases of the states |jm) within a 
given j manifold should be chosen such that the states 
follow the standard sign convention, i.e., J±|jm) = 
VjU + 1) — m{m zb l)\jm zb 1). For the case with time- 
reversal symmetry, the transformation matrix that diago- 


The ^-signature is the quantum number (ibl) that corresponds 
to the discrete symmetry transformation . 

We note that the reduced matrix element 0110^^^117) is negative 
so that (JJIQq ’\jj) < 0 is the lowest eigenvalue. Indeed, in 
an orbital with good quantum numbers nlj, {nljj\QQ ^\nljj) = 
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nalizes Qg ' is real. Thus each orbital \ jm) is determined 
up to an overall sign. 

These signs can be determined by calculating the ma- 
trix elements of Qj.{ and requiring their sign to be con¬ 
sistent with the signs determined by the Wigner-Eckart 
theorem 

= 1:Z^^0-||q(2)||j') (17) 

where < 0. Starting from the state \m\ = j 

and using the ladder operators we can successively 
descend to states with |m|=j—l,...,l/2. At each step 
we determine the correct sign of the state to be consis¬ 
tent with the sign of the corresponding Clebsch-Gordan 
coefficient in m- We note that since only positive 0 - 
signature orbitals are used, the sign of m alternates in 
the above process and therefore we need to calculate ma- 
trix elements of the type {jm\Q\^{\jfh!) (where to is the 
time-reversed state of to). 

In practice, we calculate the matrix elements 
{ja\Q^^{\jb) in the original basis \ja) obtained with ev8 

(i.e., before diagonalizing Qg '^) and then transform both 
|_)a) and \jb) to obtain the desired matrix elements 
{jm\Q^^l\jrh') (note that the same transformation ma¬ 
trix applies to both \ja) Ijto) and \jb) |jfh)). 

The ev8 code economizes on wave function storage by 
using only the amplitudes associated with the first oc¬ 
tant, x,y,z > 0. The symmetry properties of the wave 
functions under the x = Q, y = 0 and z = 0 plane re¬ 
flections should then be taken onto account in the eval¬ 
uation of the quadrupole matrix elements; see the Ap¬ 
pendix for details. In the code and its output, the orbital 
wave functions are represented on a lattice in coordinate 
space. Once we have determined the standard spherical 
orbitals (jjjm and the deformed orbitals (j)k, we can cal¬ 
culate the overlap matrix Ukjm = by taking the 

scalar product of the respective orbitals on the lattice. 
The spherical-basis matrices of the (deformed) ground- 
state density matrix p and the single-particle Hamilto¬ 
nian h are calculated from © and 0, respectively. 

V. 2®SI IN THE SCMF 

We now carry out the mapping for the SCMF of 
^®Si following the procedures described in Sect ions |TTI and 
ESI We use the SCMF code ev8, which takes as input 
a Skyrme parameterization of the mean-held energy and 
a possible pairing interaction. In this work we use the 
SLy4 parameter set 0. The hnal energies in this work 
are calculated in the HF approximation, with no pairing 
contribution.^^ 


As discussed earlier, we may include a pairing interaction when 
searching for the spherical solution but the energy is extrapolated 



sence of pairing with the SLy4 interaction. 


The energy landscape of ^®Si calculated using ev8 with 
a constraining mass quadrupole held is shown in Fig. 
Qualitatively, the landscape is very similar to that found 
for the USD interaction in Fig. ^ The ground state at 
the minimum of the curve has an oblate deformation with 
a value (Q) = —89.2 fm^. This is larger than the value 
found for the USD interaction, which is to be expected 
because the USD is a dehned in a restricted space. There 
are two other stationary points in the landscape, namely 
the spherical saddle point and and a shallow prolate min¬ 
imum at (Q) « 30 fm^. The deformation energy of the 
oblate solution = E^{ — Egph ~ 1.71 MeV is lower 
than the deformation energy found for the USD interac¬ 
tion. 

A standard spherical basis \jm) is constructed as de¬ 
scribed in Section csi The reduced matrix elements of 
the mass quadrupole operator Q in the spherical basis 
are then extracted from 

UWQWf) = I ^(2 0 j'TO|jTO)(jTO|Q|/TO) . (18) 

It is interesting to compare the reduced matrix ele¬ 
ments m to those of the simple spherical Woods-Saxon 
(WS) model, taking the parameters of the Woods-Saxon 
plus spin-orbit potential from Bohr and Mottelson 0. 
Fig.Elshows the ratio |(j||Q||j')sLy 4 /(j||Q||/)ws| of the 
quadrupole reduced matrix elements in the Sly4 SCMF 
and in the WS model versus |(j||<3||/)ws|- We only 
show results for the valence sd shell composed of the 
0 t? 5 / 2 ,lsi /2 and 0 ^ 3/2 orbitals. The solid circles corre¬ 
spond to diagonal elements {j = j') and open triangles 
to off-diagonal elements. The values of this ratio vary in 
the range ^ 0.95 — 1.1 (dotted lines in Fig. EJ- Thus, 


to the limit of no pairing. 












quadrupole renormalization effects in shell-model calcu¬ 
lations are more influenced by the truncation than by the 
particular treatment of the spherical mean field. Let us 
next examine those truncation effects. 



^ 9 " 0.8 
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FIG. 5: Ratio |(j| [Q||i')sLy4/0'||(3llj')ws| of the mass 

quadrupole reduced matrix elements in ^®Si versus 
KillQIbOwsI- Solid circles: diagonal elements, open 
triangles: off-diagonal elements. 


A. Truncation effects and operator rescaling 

As mentioned earlier, the transformation from the de¬ 
formed mean-field orbitals to the spherical basis is not 
unitary owing to truncation effects. The severity of the 
truncation can be assessed in several ways. The first is 
the trace of the truncated single-particle density matrix, 
Eq. ® , which should equal the number of particles in the 
truncated theory. With the truncation excluding orbitals 
beyond the N,Z = 20 shell closures, we find the num¬ 
ber of particles to be 27.75 compared to 28 for the full 
density matrix. When truncated to the valence sd shell 
(i.e., without including the Osi/ 2 , ^Pz /2 and Opi /2 core or¬ 
bitals), the trace of the density matrix is 11.89 compared 
to 12 valence particles in an sd shell-model theory. A 
more severe test of the truncated single-particle density 
matrix is Eq. ®, requiring that its eigenvalues be zero 
or one. Thus we diagonalize the truncated proton and 
neutron density matrices and examine the eigenvalues. 
For the truncation to the 12-1-12 orbits of the sd shell, we 
find that all eigenvalues are within 1.5% of the required 
values. We conclude that the projected density 
satisfies to a good accuracy the HE condition for a Slater 
determinant. 

While the particle number is rather robust to trunca¬ 
tion, the same is not true of the matrix elements of the 
quadrupole operator. To see how the effect of trunca¬ 
tion evolves, we impose an upper cutoff energy on the 
spherical orbitals and examine how the quadrupole mo¬ 
ment changes with this cutoff. This is shown in Fig. 



FIG. 6: Quadrupole moment of the SCMF ground state of 
^®Si using the SLy4 Skyrme energy functional and different 
truncations in the spherical basis. The horizontal axis is the 
maximum single-particle energy in the truncated space. Each 
solid circle corresponds to a spherical single-particle state of 
protons or neutrons. The untruncated moment of -89.2 fm^ 
is shown by the arrow on the right. 

describing |(Q)| as a function of the maximal single¬ 
particle energy used to define the spherical model space. 
With a cutoff of -1-15 MeV unbound orbital energy, the 
quadrupole moment is reduced from its full value of -89 
fm^ by about 7%. Excluding all positive energy orbitals, 
the moment is -69 fm^, a reduction of ~ 25%. The trun¬ 
cation that excludes all orbitals above the N, Z = 20 
shell closure yields {Q) = —54.4 fm^. Finally, we find a 
slightly lower value of {Q) = —52.4 fm^ when truncat¬ 
ing to the valence sd shell. Thus, within a shell-model 
theory in the truncated valence sd shell, the quadrupole 
operator should be rescaled by 89.4/52.4 1.71. 

B. Mapped Hamiltonian 

To obtain the mapped Hamiltonian, we first calculate 
the transformation matrix U between the deformed and 
spherical bases, using the output orbital wave functions 
from ev8. We then find the deformed single-particle 
density matrix Eq. and single-particle Hamiltonian 
Eq. 0 in the spherical basis. Next, we perform the mul¬ 
tipole decomposition of the single-particle Hamiltonian 
h = h(0) + h(2) using m and the Wigner-Eckart theo¬ 
rem (to convert the reduced matrix elements to matrix 
elements) to extract effective spherical shell orbital en¬ 
ergies e® (see Eq. O) and and an effective quadrupole 
field . From this point on the procedure is identical 
to what we did for the USD Hamiltonian. We construct 
the effective CISM Hamiltonian cu with a coupling con¬ 
stant g. A crude estimate for g in the Hartree approxi¬ 
mation is 5 « — 1/(/iq) = 0.046 MeV/fm^. To determine 
g more precisely we match the deformation energy of the 
effective CISM Hamiltonian IIH to its value 1.71 MeV 
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found in the SCMF theory (of SLy4). The deformation 
energy of the CISM Hamiltonian for a given value of g is 
calculated from the spherical and deformed (oblate) so¬ 
lutions of CD in the HF approximation. This procedure 
gives us 5 = 0.0594 MeV/fm^. 

The mapped Hamiltonian can be tested by solving it 
in the HF approximation and comparing with the SCMF. 
The deformation energy is the same by construction. We 
obtain a quadrupole moment of —55 fm^, which is ~ 5% 
larger than the SCMF quadrupole moment (projected on 
the sd shell). We also find the spherical occupations of 
the deformed HF ground state of the CISM Hamiltonian 
to be rather close to the spherical occupations of the 
deformed SCMF minimum (see Table IV below). 





FIG. 7: Left: versus ^ in ^^Si. Right: versus 

C The density matrix corresponds to the single-particle 
(Hartree) Hamiltonian hefi = -|- . 


C. Hartree approximation 

In the previous section, we solved the CISM Hamilto¬ 
nian in the HF approximation, i.e. including exchange. 
In this section we briefly examine the simpler Hartree 
approximation, because it allows some insight into the 
parametric sensitivities. Neglecting the exchange inter¬ 
action, the dependence of the deformation energy E^ef 
on the strength parameter g can be conveniently deter¬ 
mined in two steps. As a first step, the single-particle 
Hamiltonian 

heff (19) 

is solved as function of a coupling parameter say in 
the range 0 < ^ < 2. One needs to keep the expectation 
values of and as a function of calculated with 
the ground-state density matrix of (Hllll . For the sec¬ 
ond step, one determines the ground-state energy E{g) 
by finding the minimum in 

E{g) . (20) 

The minimization in Eq. iini) with respect to ^ is equiv- 
alent to the self-consistent condition g{hQ ')j = The 
deformation energy is then given by 

E^^fig) = EiO) - E{g) (21) 

Fig. [ 7 | shows (^^°^)^ (left panel) and (right 

panel) versus The value of ^ that minimizes the 
Hartree energy is determined through the competition 
between the two terms in Eq. CT . Performing next the 
minimization to get E{g), we find the deformation en¬ 
ergy curve shown in the left panel of Fig. |51 The panel 
on the right in Fig. 0 shows the quadrupole moment of 
the solution as a function of g. The discontinuity of {Q) 
at 5 « 0.037 MeV“^ describes a first order spherical to 
deformed shape transition. Note that the deformation 
energy E^etig) in the deformed phase is a steep function 
of g; thus g is well constrained by the fitting procedure. 
The solid triangles in Fig. |S1 indicate the SLy4 SCMF 




FIG. 8: Left: deformation energy versus coupling parameter 
g in ^®Si when solving the effective Hamiltonian o in the 
Hartree approximation. Right: quadrupole moment (Q) ver¬ 
sus g in the Hartree approximation using the same effective 
Hamiltonian. The solid triangles describe the SCMF solution 
of the SLy4 energy functional (note that the triangle in the 
left panel is a fit). 

point (in which the value oi g = 0.0466 MeV“^ is fitted 
to give the correct deformation energy). 

The Hartree deformation energy E^ef versus {Q) in the 
truncated sd space is plotted in Fig. El The solid triangle 
is the SLy4 SCMF point. It lies very close to the curve ex¬ 
tracted from the effective Hamiltonian, showing that the 
effective theory closely reproduces the quadrupole mo¬ 
ment in the truncated space. 

D. Correlation energy 

With the mapped shell-model Hamiltonian in hand, 
we are now in a position to compute the correlation en¬ 
ergy as defined in Eq. 113). The HF energy of OD is 
found to be E^^f = —147.44 MeV. The ground-state en¬ 
ergy Egs is calculated using the shell-model code oxbash, 
putting in the orbital single-particle energies and inter¬ 
action matrix elements from the mapped Hamiltonian 
(tTTll . We find Egg = —153.93 MeV. Anticipating calcu¬ 
lations in much larger shell-model spaces, we have also 
computed Egg with the SMMC. Taking /3 = 6 MeV“^ 
and a time slice of A/3 = 1/32 MeV“^ we find (us- 
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<Q> (fm^) 

FIG. 9: Deformation energy versus quadrupole moment of 
^®Si using the Hartree approximation for the effective Hamil¬ 
tonian GB in the sd-shell. The triangle shows the SCMF 
solution for the SLy4 energy functional where the quadrupole 
moment is calculated in the truncated space {sd shell). 


ing a stabilization method for the one-body propagation) 
Egs = —153.87 ± 0.07 MeV, in agreement with oxbash. 
For a time slice of A/3 =1/64 MeV“^ the SMMC result 
is Egs = -153.94 ± 0.07 MeV. 

The quadrupolar correlation energy we find for ^®Si is 
then Ecott ~ —4.55 MeV. This is quite comparable to the 
value we obtained from a similar reduction of the USD 
Hamiltonian (4.3 MeV). This gives one some conhdence 


that the quadrupole correlation energy has a magnitude 
of about 4-5 MeV in ^®Si. It is interesting to compare 
this with the quadrupole correlation energy obtained by 
a completely different approach, the generator coordinate 
method (GCM) M U sing that method in a global study 
of correlation energies [l^, the authors of Refs. [iM 
obtained a correlation energy of 4.9 MeV for ^®Si, which 
is somewhat larger. However, these calculations included 
a pairing interaction, which we have not treated here. In 
particular, the J = 0 wave functions used as the GCM 
basis are projected from deformed states that are con¬ 
structed in the presence of a delta pairing force. 


VI. OTHER sd SHELL NUCLEI IN THE SCMF 

We have repeated the mapping discussed in Section M 
for the SCMF of all other deformed N = Z nuclei in the 
sd shell, i.e., for ^°Ne, ^"^Mg and ^®Ar (^^S is found to 
be spherical, just as for the USD interaction). The results 
are summarized in Tables H, HI and IV (including results 
for 28Si). 

Truncation effects on the deformed SCMF density ma¬ 
trix are summarized in Table H. The eigenvalues and 
trace of the truncated density matrix (projected 

on the valence sd shell) are listed separately for protons 
and neutrons (there are 6 positive z-signature eigenvalues 
for each type of nucleon). We see that the HF condition 
= p holds to a rather good accuracy for all four nuclei 
(all eigenvalues are very close to either I or 0). The con¬ 
dition that the trace is given by the number of valence 
particles is also satisfied quite well. 


TABLE II: Eigenvalues and trace of the truncated density matrix in the SCMF theory of the SLy4 Skyrme energy 

functional. Shown are results for deformed 4N nuclei in the sd shell (separately for protons and neutrons). 


Nucleus 


eigenvalues of Pp^ifP 

tr (PPmfP) 

2°Ne 

protons 

0.974 

0.003 

9 TO"® 

0 

-3 TO”® 

-6T0"® 

1.953 


neutrons 

0.975 

0.003 

4-10"® 

0 

-4 -10“® 

-6 TO”® 

1.955 

^‘‘Mg 

protons 

0.977 

0.975 

0.006 

5-10~® 

-2T0"® 

-8T0"® 

3.915 


neutrons 

0.977 

0.975 

0.006 

7T0“® 

-2T0"® 

-610"® 

3.916 

28Si 

protons 

0.991 

0.990 

0.987 

0.004 

7-10"® 

-6T0"® 

5.944 


neutrons 

0.991 

0.991 

0.987 

0.004 

2 TO"® 

-8T0-® 

5.945 

36 Ar 

protons 

0.999 

0.998 

0.997 

0.995 

0.993 

0.001 

9.966 


neutrons 

0.999 

0.998 

0.997 

0.995 

0.994 

0.001 

9.966 


Table HI compares results for the SCMF with the SLy4 
Skyrme force and for the HF of the mapped CISM Hamil¬ 
tonian m with g determined by matching the deforma¬ 
tion energy of both theories. Shown are the spherical 
and mean-field total energies, deformation energy and 


the quadrupole moment {Q) in the truncated sd shell- 
model space. We observe that the quadrupole moment 
of the effective theory deviates by at most 5% from its 
value in the SCMF (for ^Sgi). 

We also show in Table HI the total SCMF quadrupole 
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TABLE III: SCMF with the SLy4 Skyrme energy functional and the mapped Hamiltonian Eq. llTTl for deformed 4A^ sd shell 
nuclei. All energies are in MeV, (Q) is in fm^ and g is in MeV/fm'^ 


Nucleus 

Interaction g 

-E^sph -^mf -^def (Q) (Q)mf -E^corr 


SLy4 

0.077 

-151.83 -157.43 5.6 49.0 84.0 

-30.17 -35.78 fit 48.7 2.39 

'"Mg 

SLy4 

0.0394 

-187.20 -195.92 8.73 59.6 111.9 

-79.29 -88.03 fit 60.6 3.50 

2«Si 

SLy4 

0.0594 

-233.19 -234.90 1.71 -52.4 -89.2 

-147.67 -149.38 fit -55.0 4.55 

36 Ar 

SLy4 

0.117 

-303.16 -305.43 2.27 -43.4 -74.7 

-274.53 -276.79 fit -42.7 1.42 


TABLE IV: Spherical occupations (hj) for deformed AN nuclei in the sd shell. The SCMF occupations (computed for the 
SLy4 Skyrme energy functional) are compared with the occupations for the mapped Hamiltonian Eq. Hill both in the HE 
approximation and in shell-model calculations. 


Nucleus 

Interaction 

P 0^5/2 

P lSl/2 

P 0^3/2 

n 0^5/2 

n lsi /2 

n 0d3/2 


SLy4 

1.481 

0.376 

0.096 

1.486 

0.374 

0.095 

2°Ne 

(HF) 

1.565 

0.372 

0.064 

1.569 

0.366 

0.065 


(CISM) 

1.531 

0.403 

0.066 

1.539 

0.396 

0.065 


SLy4 

3.290 

0.411 

0.215 

3.293 

0.408 

0.215 

'"Mg 

(HF) 

3.398 

0.346 

0.256 

3.399 

0.344 

0.257 


(CISM) 

3.356 

0.329 

0.315 

3.360 

0.327 

0.313 


SLy4 

4.956 

0.723 

0.265 

4.965 

0.717 

0.263 

to 

00 

(HF) 

4.929 

0.712 

0.359 

4.939 

0.702 

0.359 


(CISM) 

5.011 

0.633 

0.356 

5.010 

0.633 

0.357 


SLy4 

5.919 

1.751 

2.297 

5.921 

1.758 

2.287 

36 Ar 

/i(2) (HF) 

5.979 

1.694 

2.327 

5.980 

1.702 

2.3185 


(CISM) 

5.971 

1.662 

2.367 

5.971 

1.668 

2.361 


moment ((3)mf calculated in the deformed HF ground- 
state. Comparing {Q)mi with (Q) in the trun¬ 
cated sd shell, we estimate that in the CISM the¬ 
ory the quadrupole operator should be rescaled by « 
1.71,1.88,1.70 and 1.71 for 20Ne,24Mg,28Si and ^SAr, 
respectively. 


The correlation energies (calculated from Eq. 11511 with 
Egs found using oxbash) are shown in the last column of 
Table III. The values we find are quite reasonable. They 
are below the values cited in but we have not yet 
included pairing in our effective theory. 


Finally, in Table IV we compare the occupations of the 
spherical valence orbitals of the SLy4 SCMF theory with 
similar occupations obtained with the mapped Hamilto¬ 
nian CD- The latter are calculated in the HF approx¬ 
imation as well as in the shell-model theory from the 
expectation values of h, = in the correlated 

CISM ground state. 


VII. THE Q Q INTERACTION 

In nuclear physics there is a long history of model¬ 
ing the effective interaction as a quadrupole-quadrupole 
interaction of the form Q ■ Q, including Elliott’s SU(3) 
model 13 and the pairing-plus-quadrupole model of 
Kisslinger and Sorensen [3- Since it is possible in some 
contexts to exploit the properties of the quadrupole op¬ 
erator, it is of interest to see how well a quadrupole- 
quadrupole interaction performs in our context. Thus, 
we apply the Hamiltonian 

H = -^gQ-Q. (22) 

in the truncated shell-model space with the matrix ele¬ 
ments of and the mass quadrupole operator Q deter¬ 
mined by the appropriate mean-field theory. We follow 
the same procedure we used before to determine the cou¬ 
pling constant g. 

We have determined such an effective Q ■ Q 
interaction for the deformed AN sd shell nuclei 
^°Ne, ^^Mg, ^®Si and ^®Ar , using the HF approximation 
of the USD interaction as the SCMF theory. We observe 
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that the quadrupole moments calculated with both ef¬ 
fective interactions m and (P|) are quite close to the 
USD results. This can be seen in Fig. EH where (Q) 
(in fm^) is plotted versus mass number A. Results for 
the Q ■ Q interaction (solid triangles) and the interaction 
with (solid squares) are compared with the USD re¬ 
sults (open circles). The effective Hamiltonians (Hill and 
(P|l also yield similar occupations {fij) (see Fig. 121) and 
correlations energies (see Fig. E]). 



FIG. 10: Quadrupole moment in the SCMF method versus 
mass number A for AN sd-shell nuclei. Results are shown 
for two effective CISM Hamiltonians: the Q ■ Q interaction 
II22II (solid triangles) and the interaction Ullll (solid 

squares). They are compared with results for the USD Hamil¬ 
tonian (open circles). 


VIII. PROSPECTS AND CONCLUSIONS 


We are encouraged by the results we found for con¬ 
structing a isoscalar quadrupolar interaction from the 
SCMF in the sd-shell, but much remains to be done to 
demonstrate that the method will be useful in deriving a 
global theory of correlation energies. We list some of the 
major tasks that should be addressed: 


placed by 


hL = 


'^ac y l^M\ ) "'M-ac\l^bd"'-M-bdf^db) 


+5Sm( J2bd^^M-,ad^-M-,bcPdb ^Qac , (23) 


where is the self-consistent density matrix in 
the spherical basis (here a = (jm)) and Q is the 
axial quadrupole operator. A practical problem is 
how large should we take A to be. One possibility is 
to use the value of A that gives the deformed J = 0 
minimum. 


• The real challenge of a global theory is the heavy, 
open-shell nuclei, and we should ask whether any 
severe problems will arise when we attempt to map 
the SCMF onto the CISM Hamiltonian. As a con¬ 
crete example, consider the nucleus ^®^Dy. It is 
strongly deformed, having an intrinsic quadrupole 
moment around 18 bn. The SCMF using the SLyd 
interaction reproduces this value very well |l8| . 
The deformation energy when pairing is not in¬ 
cluded is quite large (about 20 MeV). The J = 
0 projection and quadrupole fluctuations (in the 
presence of a delta pairing force) contribute an ad¬ 
ditional correlation energy of « 3.4 MeV [T^. 

To use an CISM Hamiltonian to calculate defor¬ 
mation and correlation energies, one would need 
to truncate to a basis of the full shells, i.e., N = 
82 — 126 for neutrons and Z = 50 — 82 for pro¬ 
tons plus possibly some intruder states. Thus, a 
first test would be to see how well the consistency 
conditions Eqs. and are satisfied in such a 
truncated space. 

The next question that arises is the practicality of 
various implementations of the CISM Hamiltonian 
for solving in such large spaces. There are several 
approaches that have promise, including the shell 
model Monte Carlo (SMMC) method 0, 1^ - 
the Monte Carlo shell model and direct ma¬ 
trix methods such as those used by the Strasbourg 
group 1^ . One could also consider more approxi¬ 
mate theories of the correlation energy, such as the 

rpaH. 


• So far, we have carried out the mapping only for 
nuclei that have deformed HF ground states. If the 
lowest-energy mean-field solution is spherical, we 
can use a constraining quadrupole field XQ to de¬ 
form the HF ground state and use this deformed 
solution to construct the effective CISM Hamilto¬ 
nian. In particular, for a given strength A of the 
field, we can compare the deformation energy of the 
constrained SCMF theory with the deformation en¬ 
ergy of the effective Hamiltonian in the constrained 
HF approximation. The Hartree-Fock Hamiltonian 
da of the effective CISM Hamiltonian is now re- 


• The mapping should be extended to include all 
other parts of the interaction that we believe are 
well-determined in the SCMF, in particular the in¬ 
teraction in the isovector monopole and quadrupole 
densities. This should be straightforward to carry 
out, requiring only that nuclei with N ^ Z he con¬ 
sidered in a more general parameterization of the 
interaction. 


Nuclei with spherical HF mini mum acquire deformation after 
applying a J = 0 projection M- 
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• The pairing interaction is very important in nuclear 
structure but is notoriously difficult to specify in 
detail, i.e., its density dependence and its spatial 
range (or equivalently the orbital space energy cut¬ 
off). Most SCMF energy functionals include a pair¬ 
ing interaction, but only an overall strength can be 
reliably determined in phenomenological theories. 
Thus, absent an ab initio theory of the pairing in¬ 
teraction, we cannot expect to have a reliable de¬ 
composition of the pairing fields in the mean-field 
Hamiltonian. Still, it should be possible to deter¬ 
mine a strength parameter controlling the overall 
magnitude of pairing effects. Note that in the pres¬ 
ence of pairing, the single-particle density matrices 
of the SCMF theory have the BCS form. Thus, to 
determine the parameters in the mapped Hamilto¬ 
nian, one would use the HFB or the HF-BCS ap¬ 
proximations for the self-consistent mean-field the¬ 
ory. 

• An unsatisfactory aspect of our mapping is the lack 
of unitary of the transformation from the deformed 
basis to the truncated spherical basis. It would be 
useful to investigate mapping formalisms that yield 
density matrices that satisfy exactly the criteria we 
have discussed, i.e., Eqs. and ©■ 

Despite this long list of tasks to be done, we believe 
that the method discussed here can be developed to pro¬ 
duce a well-justified global theory of shell-model correla¬ 
tion energies, starting from a global SCMF theory (e.g., 
a Skyrme energy functional). Particularly encouraging 
is the finding that quadrupole correlations account for 
more than half the correlation energy in sd shell nuclei. 
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Appendix: wave functions and matrix elements in 
the 2 -signature representation 

A deformed single-particle wave function in the density 
functional theory is a two-component complex spinor '0 = 


( alternatively a four-component real spinor ^|J = 

*3 with tpi = Re 'ipj- ; D2 = Im ; -03 = Re '0_ and 

04 = Im0_. The SGMF orbitals are chosen to have good 
parity tt and good 2-signature ^ . Such states have well- 
defined symmetry properties under each the three plane 
reflections x —x; y —y; z —z (see Table I in 
Ref. 0 ). Therefore it is sufficient to store their values 
in one octant (e.g., x,y,z > 0). 

The scalar product of two spinors 0 and 0 is a complex 
number with Re (0^0) = I™ (0^0) = 0i02~ 

0201 + 0304 - 0403- 

The find a standard basis of spherical single-particle 
orbitals \jm}, it is necessary to calculate matrix elements 
of the quadrupole operators (see Section Hg. We 
denote the j -I- 1/2 orbitals with good j and positive 2- 
signature by \ja). The quadrupole matrix elements are 
obviously real and given by 


{ja\Qo^ \jb) =8 dxdydz {2z‘^-x‘^-y‘^) Re (0'l'0) , 

J a:,y,2:>0 

(24) 

where 0 and 0 are the spinors describing \ja) and \jb), 
respectively. 

To calculate {ja\Q^^{\jb) we use the symmetry proper¬ 
ties of the orbitals under plane reflections and find 


{ja\xz\jb) =8 dx dy dz xzKe (f>) 

J x,y,z>0 

{ja\yz\jb) = 8i j dxdy dzyzluviiP^) , (25) 
Since the time-reversed spinor of 0 is given by 0 = 



Re(0'f0) = 0103 - 0301 +0402 - 0204 
Im (0'f0) = 0302 - 0203 + 0401 - 0104 ■ (26) 

The matrix elements of (§±1 ^'I'e calculated from 12511 
using Q±i = +V6l]0a:i2i ± iyiZi) (they are all real). 
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